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Abstract 

Within the mode-coupling theory (MCT) for the evolution of structural relax- 
ation in glass-forming liquids, correlation functions and susceptibility spectra 
are calculated characterizing the rotational dynamics of a top-down symmet- 
ric dumbbell molecule, consisting of two fused hard spheres immersed in a 
hard-sphere system. It is found that for sufficiently large dumbbell elonga- 
tions, the dynamics of the probe molecule follows the same universal glass- 
transition scenario as known from the MCT results of simple liquids. The 
a-relaxation process of the angular- index-j=l response is stronger, slower 
and less stretched than the one for j=2, in qualitative agreement with re- 
sults found by dielectric-loss and depolarized-light-scattering spectroscopy for 
some supercooled liquids. For sufficiently small elongations, the reorienta- 
tional relaxation occurs via large-angle flips, and the standard scenario for 
the glass-transition dynamics is modified for odd-j responses due to precur- 
sor phenomena of a nearby type-A MCT transition. In this case, a major 
part of the relaxation outside the transient regime is described qualitatively 
by the /3-relaxation scaling laws, while the a-relaxation scaling law is strongly 
disturbed. 
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I. INTRODUCTION 



During the past ten years, the evolution of structural relaxation in glass-forming liquids 
has been intensively studied using neutron-scattering spectroscopy, various light-scattering 
techniques, dielectric-loss spectroscopy, and molecular-dynamics simulation. Results of this 
work have also been used to test the mode-coupling theory (MCT), which interpretes the 
structural relaxation as precursor of the glass-transition. Originally, the MCT was proposed 
as an approximation approach for the cage effect in liquids In its simplest version, 

the MCT equations of motion describe an ideal liquid-to-glass transition, i.e. a bifurcation 
from ergodic to nonergodic dynamics, if control parameters like temperature T or packing 
fraction ip cross critical values T c or ip c , respectively. This bifurcation is connected with the 
evolution of a two-step relaxation scenario entirely determined by the regularly changing 
equilibrium structure. Two divergent time scales appear, closely connected to two power- 
law decay processes. A detailed description of these results can be found in Ref. |J and 
references therein. Comparisons of the theoretical results for simple model systems with 
experiments done on colloids j|,|5|], and with computer-simulation studies |||7j demonstrate 
the validity of the microscopic MCT approach. For the solutions of the MCT equations, a 
variety of results has been derived by asymptotic expansions, using as a small parameter 
the distance from the critical point, e = (ip — ip c )/ip c , or e = (T c — T)/T c , respectively. The 
leading-order results of this expansion establish universality features of the MCT dynamics. 
Assessments of the theory have been reached by comparing spectra in the GHz regime or 
relaxation curves within the pico-second window with the universal results. The outcome of 
this work, which is reviewed in Ref. || , leads to the conclusion that MCT properly describes 
some essential features of structural relaxation even for some complicated molecular liquids. 

The MCT for simple systems has been extended recently to liquids of nonspherical 
molecules PHTT|. But so far, only the bifurcation equation for the so-called nonergodic- 
ity parameters resulting within the new theory could be solved. Comparing these results 



with the findings of molecular-dynamics simulations for a liquid of linear molecules [12,13 
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and for water |TT[ indicates that the MCT for molecular liquids is promising. It was also 
predicted that there can be two states of nonergodic motion for nonspherical molecules. 
These states are connected by a type-A transition if the molecules exhibit a top-down sym- 
metry P, |l0| , [i4|j . Such transitions are generic possibilities in MCT, provided there is some 



symmetry in the problem rendering certain mode-coupling coefficients zero [|15|| . At a type-A 
transition, the nonergodicity parameters change continuously, whereas at the conventional 
MCT transition, referred to in this context as a type-B transition, a discontinuity occurs 
Of. 

In this paper, correlation functions and susceptibility spectra shall be discussed, which 
deal with the glassy dynamics of the orientational degrees of freedom of nonspherical 
molecules. The results are obtained as solutions of the equations of motion derived pre- 
viously for the dynamics of a linear probe molecule immersed in a simple liquid. A 



top-down symmetric dumbbell of two fused hard spheres will be considered as the molecule, 
and as the solute, a hard-sphere system is chosen. This model deals with the simplest 
problem concerning glassy rotational dynamics, namely the influence of the cages formed 
by the neighbors of the molecule on the molecule's reorient at ional motion as it is caused by 
steric hindrance. The dynamics will be exemplified for two molecule with a large 

elongation and a molecule with a small elongation. 

It will be shown that large elongations lead to strong coupling of the rotational degrees 
of freedom to the density fluctuations of the solute, such that the glassy dynamics of the 
latter enforces the validity of all the universal MCT laws for the solvent. Moreover, the 
corrections to the leading-order-asymptotic laws show the same qualitative trends as studied 



for simple liquids |17|,|18[ . A motivation of the present study is the explanation of three 
general properties of the a-relaxation in molecular liquids, which are exhibited in Fig. [I]. 
In this figure, experimental susceptibility spectra for the van-der-Waals liquid propylene 
carbonate (PC) are reproduced for four temperatures. One set of data deals with the 
response for angular- momentum index it was obtained by dielectric-loss spectroscopy 



19fl . The other set was measured by depolarized-light-scattering spectroscopy |2(| and deals 



with the (.7=2) reorient at ional dynamics. The data show for T = 293K and T = 295K a- 
relaxation peaks at 4GHz (j = l) and 10GHz (j=2) respectively. These temperatures exceed 
the melting temperature T m = 218K of PC by more than 70K. Lowering T to 200K, the 
a-peaks of the spectra are shifted down by about two orders of magnitude. The shape of the 
a-peak is temperature independent, and the ratio of the a-pro cess-time scales, characterizing 
the a-peak-maximum positions for the two values of j, is also T-independent. These are 
two features which MCT predicts to be universal. The first nonuniversal feature to be 
understood is, that the a-peak intensity, taken relative to that of the band of microscopic 
excitations at around ITHz, is larger for the (j=l) response than for the (j=2) case: the 
former exceeds the latter by about a factor of 2.7. Second, the (j=l) response is slower than 
the response for j=2: the ratio of the a-peak positions is about 2.5. Third, the a-peak of 
the (j=l) response is less stretched than the peak for j=2, i.e. the halfwidth of the (j=l) 
peak is smaller than that of the (j=2) peak. If one describes these peaks by the spectra of 
the Kohlrausch law, $(£) oc e"^ T ^ , the stretching exponent (3 for j=l, f3j=i ~ 0.9 JOJ, is 



larger than the one for j=2, (3j =2 ~ 0.8 The same three a-peak features are noticed, if 
one compares the depolarized-light-scattering spectra of glycerol jnj] with the corresponding 
dielectric-loss spectra |22| . A fourth general feature to be explained is the large ratio of the 



a-relaxation-time scale found by depolarized-light-scattering spectroscopy and the one found 
for the longitudinal elastic modulus by Brillouin-scattering spectroscopy. For Salol, a ratio 



of about 10 was reported ||23|| , while for PC, a factor of about 5 was found [^0 . 

The small elongation of concern in this paper is chosen so that it exceeds the critical value 
for the above-mentioned type- A transition by about 10%. The theory for the corrections to 
the leading-order asymptotic laws |H],E01 i m phes that these diverge at a type-A transition. 
Therefore, the range of validity of the universality features of the standard MCT bifurcation 
shrinks upon approaching the type-A transition. It will be shown that in our example the 
standard results are not exhibited any more for reasonable choices of the distance parameter 
e. In particular, it is impossible to identify a two-step scenario for the odd-/ correlators, nor 
is a-relaxation scaling observed. 



The paper is organized as follows. In Sec. 0, the model is defined, and the MCT 
equations are noted. After an overview of the general scenario for the evolution of the 
glassy relaxation of the reorientational correlators (Sec. IIIA), the differences between the 
relaxation patterns for the (j=l) and (j=2) response are described for strong (Sec. IIIB) 
and weak (Sec. IIIC) steric hindrance. In Sec. HID it is demonstrated how the /3-relaxation 
is described by the first scaling law, and in Sec. HIE it is discussed how the a-relaxation- 
scaling-law description emerges. The concluding Sec. IV summarizes the results. 



II. THE MODEL SYSTEM 



A. The solvent 

A system of N spherical particles shall be considered as the solvent. The basic variables 
describing the structure are the density fluctuations for the wave vector q: g$ = ^2 K exp(iq • 
r K )/v / A. Here r* K , k = 1,2,.., A, labels the centers of the particles. The structure factor 
<3q — (\Qg\ 2 ) provides the simplest information on the equilibrium distribution of the particles; 
here, ( ) denotes canonical averaging. Because of rotational symmetry, S q depends on the 
wave- vector modulus q = \q\ only. The structure factor can be expressed through the direct 
correlation function c q via the Ornstein-Zernicke equation S q = 1/(1 — pc q ); where p denotes 
the particle density fl2"l| . The simplest quantities, characterizing the structural dynamics in a 
statistical manner, are the normalized auto-correlation functions for the density fluctuations, 
called the density correlators $ q (t) = (Qq{t)* Qq*) / S q . The evolution with increasing time t is 
given by the canonical equations of motion. We will also need Fourier-Laplace transforms 
for complex frequency z, Imz > 0, ® q (z), using the convention: F(z) = i J °° exp(i;zt) F(t)dt. 
For real frequency u, one gets with z = u + iO: F(z) = F'(u) + iF"(u). The imaginary part 
F"(uj) is called the fluctuation spectrum, and x"(uj) = ujF"(uj) is the susceptibility spectrum 



The basic version of MCT consists of two equations j2| . The first one is exact and derived 
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within the Zwanzig-Mori formalism: 

d 2 $ q (t) + n 2 q $ q (t) + n 2 q f m q (t - O d t ,® q {t') df = . (la) 

Jo 

Here, f2 9 = vq/ y/~S~ q , with t> denoting the thermal velocity. The relaxation kernel m q (t) 
is a fluctuating-force correlator. The equation has to be solved with the initial condition 
$ 9 (t) = 1 - {Vt q t) 2 /2 + C(t 3 ) [13. Equation flTa|) is equivalent to the double fraction 

= ^2 • (lb) 

z — 

z + Q? q m q {z) 

The second MCT equation is obtained by writing the kernel as a sum of a regular term and 
a contribution describing the cage effect. The latter is treated by Kawasaki's factorization 
approximation for the force correlations. It is found to be a quadratic functional of the 
density fluctuations: Y^k+p=ff V(q,k,p)^k{t) ^p(t)- For the sake of simplicity, the regular 
term shall be neglected in the following. Furthermore, the wave-vector modulus will be dis- 
cretized to M values with equal spacing h. Thus, q, k, p can be considered as labels running 
from 1 to M. As a result, the kernel is given as a quadratic mode-coupling polynomial T q 
of the M correlators <& q (t), q = 1, M: 

m q (t) = .F ff [$*(*)] = V ^p ^if) %{t). (2) 

kp 

The positive coupling coefficients V q k p are given by S q and c q |L7|]. Anticipating these equi- 
librium quantities to be known, Eqs. ( ]Ta| ) and (0) are closed. 

Equations ([!]) and (fj) exhibit a transition from liquid-state dynamics in the regime 
T > T c or ip < ip c to glass-state dynamics for T < T c or ip > (p c . In the former regime 
the density fluctuations decay to zero for long times, — > oo) = 0. The ideal-glass 
states exhibit a nontrivial long-time limit, which is called the nonergodicity parameter, 
f q = $ 9 (£ — >• oo) > 0. It is the Debye- Waller factor of the glass. At the transition, this 
long-time limit is discontinuous, and the jump is called critical nonergodicity parameter 
or plateau, f q = f q (T /* T c ,(p \ ip c ) > 0. At the critical point, the correlators decay 
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algebraically: $,(*) = fg+h q {t/t )- a + O((t/t y 2a ). The exponent a, < a < 1/2, is called 
the critical exponent, and h g > is denoted as the critical amplitude. t marks the time 
scale of the transient from the microscopic motion to the relaxation dynamics of the MCT. 
The MCT a-process is defined as the dynamics for those times, where the correlators of 
the liquid decay from the plateau fi to zero. The MCT /3-process deals with the dynamics, 
where the correlators are near the plateau, i.e. \& q (t) — f° | -C 1. The first relaxation step of 
the anomalous dynamics is given by the initial part of the /3-process; it deals with the decay 
towards the plateau for times outside the transient: (t/t ) ^> 1, & q (t) > fg. The second step 
is the a-decay in the liquid. Its initial part is identical with the final part of the /3-process, 
and it follows von Schweidler's law — f q oc —h q t b . The exponent b, < b < 1, is called 
the von-Schweidler exponent. 

In a leading-order expansion in the small parameter — fg \ one finds the universal 

results for the /3-process. There holds the factorization theorem 

%(t)-r g =hgG(t). (3) 

The dependence on time and on control parameters is given by the g-independent function 
G(t), which is called the /3-correlator. It is determined by the equation 

t 

G(t - t')G(t') dt', (4) 

to be solved with the initial condition G(t — > 0) = (t/to)~ a + 0(t a ). The number A, 
< A < 1, is referred to as the exponent parameter, a is a smooth function of the 
control parameters and is called the separation parameter. Its zero defines the critical 
point. Expanding in leading order in the distance e, one can write a = Ce , C > 0. 
From Eq. (Q), one derives the first scaling law 

G(t) =c a g±(i),e%0 t = t/t a . (5) 

Here, c a = \f\o\ denotes the amplitude scale, and t a abbreviates the first characteristic time 
scale of the MCT-transition scenario: 
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= VM 1/2a . (6a) 
The master functions g±(i) are determined by A as solutions of Eq. (M for a — ±1, respec- 



tively. They interpolate monotonously between g±(t <C 1) = t~ a and g + (t ^> 1) = — A 

or 3> 1) = —Bt b . Von Schweidler's law is obtained as the long-time limit on the scale 
tcr in the form: — fq ~ h q (t/t' a ) b . Here t^. abbreviates the second characteristic scale 

of the theory: 

C = fT^VkP , 7 = (l/2a) + (1/26). (6b) 

The a-process obeys for e — > the second scaling-law, called the superposition principle: 

$ q (t) = <f> 9 (t) , t = t/C (7) 

The control-parameter independent master function <& q (t) exhibits the initial decay = 
fq ~ h q t b + 0(t 2b ). The parameters /|, /i g and A are determined by T q from Eq. (H) for 
control parameters at the critical point. The same holds for the function The constant 

C is determined by the first Taylor coefficient in e of the deviations of T q from its value at 
the critical point. Formulas for these quantities can be found in Ref. |Tj|, where also the 
original work is cited. The theory for the leading corrections to the quoted results has been 



worked out in Ref. 17 



The calculations in this paper will be done for the hard-sphere system (HSS). The tem- 
perature does not enter the structure, but determines the time scale via the thermal velocity 
only: v 2 oc T. The relevant control parameter is the packing fraction: tp = n(pd 3 )/6, where 
d is the particle diameter. The structure factor will be calculated within the Percus-Yevick 
theory EM. The discretization will be done for M = 100 wave- vector values with step size 



hd = 0.4. For this model, all the mentioned MCT quantities have been reported in Ref. |L7 . 
In particular it was found: <£> c = 0.516, C = 1.54, 

A = 0.735, a = 0.312, b = 0.583, 7 = 2.46, B = 0.836. (8) 

The results for the glass-transition of the HSS are documented comprehensively in Refs. 



17|,[18| , albeit for a Brownian microscopic dynamics. The bifurcation scenario for the model 



with Newtonian dynamics as denned by Eqs. ([D and (Q) is demonstrated for the wave- vector 



q = 10.6/d in Ref. f25(| , where the transient time-scale was determined: to = 0.0236(d/v). 
For the presentation of our results in the following figures, the units of length and time 
will be chosen so that d — 1 and d/v = 1. The control parameters <p shall be cited by the 
logarithm x of the distance parameter e: 

(ip - <p c )/<p c = e = ±10"*. (9) 



As in the previous work the MCT equations are solved in the time domain. The 



solutions are then Laplace-transformed to get &' q {oj) + Similarly, the transformed 

kernel rnfJu) + im''(uS) is calculated from m q (t) in Eq. @. These results are used to 
compare the left-hand side of Eq. (|T5]) with the right-hand side. Thereby a verification of 
the numerical solutions is obtained. 



B. The solute 

As a model for a dilute solution of molecules we shall consider a single linear molecule 
immersed in a simple system. The position of this molecule is described by the tensor- 
density fluctuations Qj(q) = R^(e) exp(ig • f). Here, r denotes the center-of-mass position 
and e abbreviates the axis of the molecule. The RJ' are related to the spherical harmonics 
by: Rjf(e) = V yAn YJ(e). The solute-solvent equilibrium correlations are described by 
the generalized structure factors Sj(q) = (£*(<?o)f?j(<?o)), where <fo = (0,0, q). The proper 
generalization of the density correlators for simple systems are tensor-density correlators 
for the molecule, (g% (q , t)*Qj(q )). The MCT for these quantities shall be simplified by 
restricting the correlators to the diagonal elements 

(io) 

Correlation functions for wave vectors q different from go can be obtained from the specified 
ones by elementary transformations jllj. 
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The first equation of the MCT for the molecule dynamics reads [10]: 



-1 



z) = 7V2 7V2 • ( U ) 



o 2 o 2 

1 L Tq 1 c Rj 



z + Q% q (q) m T (qj/j,, z) z + fi^ m R (gj>, z) 
Here, Qx« = is the characteristic frequency for the translational motion of a tagged 



particle. = v^^/JQ + 1) is the analogue for the rotational dynamics, where j(j + 1) 

plays here and in the following a similar role as q 2 for the translational motion. The frequency 
t>R denotes the thermal velocity for the rotation. The relaxation kernels and are 
approximated along the same lines as indicated above for simple systems. They are obtained 
as a functional of the density correlators of the solvent, multiplied by the tensor-density 
correlators of the solute |10| . Let us discretize the wave vector to, say, M' values with equal 
spacing h' . Let us also restrict the angular-momentum index by some upper cutoff value l co . 
One obtains the kernels as mode-coupling polynomials 

m a (qjfx, t) = T aqjtl [$(klv, t), $ p (i)] 

= W $(klu, t) %(t) , a = R, T. (12) 

kplu 

The positive coefficients V aq jfj,(kplu), j, I = 0, 1, .., l co are given in Ref. p6f as specialization 
of the results in Ref. fL0| . They are expressed in terms of S q and Sj(q) for J = 0, 1, 2l co . 



Anticipating S q , Sj(q) and as known, Eq. (pT|) and (]T2|) are closed equations for the 

determination of the M' ■ (l co + l) 2 correlators 

The quantities of main interest for a statistical description of the rotation of the molecule 
are the reorientational correlators, defined with the Legendre polynomials Pj\ 

C^(t) = (P J (e(t)-?)}i J = 1,2,... . (13a) 

They are the long- wave-length limits of the general correlators; $(g — > j0,t) = C^\t) 
|T0| . One gets from Eq. ([TT]) the fraction representation in analogy to Eq. (|TI5|): 

C«\z) = =-5-2 . (13b) 

z + n 2 Rj mf{z) 
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Here the kernel mj(z) is the q — > limit of m-n(qjO, t). Carrying out the limit in the general 
formula for m^(qjO, t) [I(J and discretizing the wave-vector integral afterwards, one finds: 

mf(t) =rMkjiM,t),* k (t)] = V L *(klv,t)* h (t). (14) 

j klv 

The positive coupling coefficients V^ lu are listed in Ref. [p6][ . After evaluation of Q(klb>, t) and 
$fc(t) the correlators for the M' values of k, Eq. flllD yields the kernel m^(t). Fourier-Laplace 
transformation gives rr$(z) and Eq. ( |13b| ) provides C^'(z). Fourier-cosine transformation 
of the spectrum C®"(u) leads to C®(t). 

The theory shall be applied for a dumbbell consisting of two equal fused hard spheres of 
diameter d and distance (d between the centers. Thus, besides the packing fraction (p, there 
is the elongation parameter ( as the second control parameter specifying the structure. The 
structure factor Sj(q) and the corresponding pair correlation functions gj(r) are evaluated 
within the Percus-Yevick theory ||27|| . Figure || exhibits the probability distribution to find a 



solvent particle in the plane through the symmetry axis of the dumbbell. The upper panel, 
calculated for ( = 0.80, shows a pronounced quadrupolar pattern extending over several 
shells. For the small elongation £ = 0.33, the lower panel shows that anisotropy is almost 
lost from the third shell onwards. The calculations of the dynamics will be done for such 
moment of inertia that v^/d = \/2v/d. The discretization will be done with M' = 50 wave 
vectors with spacing hi = 0.8. The cutoff for the angular-momentum index is chosen as 
lio = 7 for ( = 0.80 and l\ = 5 for ( = 0.33. The equation of motion ( [ill) is transformed 



to an integro-differential equation in analogy to Eq. (|Ta|) and then solved by an algorithm 



similar to that used for the standard MCT problem [26 



III. RESULTS 



A. General features of reorientational relaxation 



Figure |3] demonstrates the transition scenario for the solute correlators for two represen- 
tative wave vectors q and three values of the angular momentum index j. The calculated 
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correlators exhibit a very weak dependence on the helicity index /i, and therefore only the 
solutions for /x = are shown. The wave vector q = 7.0 is close to the structure-factor-peak 
position, and q = 10.6 is near the first minimum of S q . The correlator for j = is the 
probability distribution of the molecule's center-of-mass position, i.e. the analogue of the 
incoherent-intermediate scattering function for simple liquids: $(g00,t) = Results 
for j=l and j=2 deal with the propagation of the dipole- and quadrupole-density fluctu- 
ations, respectively. The critical-decay curves, i.e. the solutions for ip = ip c , organize the 
bifurcation pattern. They deal with the stretched decay towards the plateaus f c (qj0). If tp 
increases above ip c , the long-time limits f(qj0) = $(gj0,t — > oo) increase above the plateau 
because the molecule gets more tightly localized in the frozen solvent. The /(gjO)-versus-g 
curves are bell-shaped, since the molecules are localized with a nearly Gaussian probability 
distribution flQfl . For ip < ip c , the correlators exhibit a long-time decay from the plateau to 



zero, and this is the a-process. The a-decay time is the larger the smaller the wave vector 
q, while the a-relaxation stretching increases with increasing q. The g-dependence of the 
relaxation features are similar as observed and explained previously for the tagged-particle 
correlator & s Jt) for simple liquids. Therefore, the following discussion shall be restricted to 
the (q = 0) limit, i.e. to the reorientational correlators C^\t). 

Figure |] exhibits representative decay curves C^\t) for the liquid state for two sepa- 
ration parameters, and Fig. |5] exhibits an extended set of susceptibility spectra x i u ) — 
u>C^'"(uj). The plateaus and critical amplitudes shall be denoted by /? and hj, respectively. 



They have been calculated from the mode-coupling functionals |10[ , and some examples are 



listed in Tab. |. These parameters specify the leading-order asymptotic results for the (3- 
relaxation process as explained in Sec. |II A| for the solvent. The factorization theorem holds 
in analogy to Eq. <g) |IJ: 

C®(t) = f? + h j G(t). (15) 

The /3-correlator G is the same function as explained in connection with Eqs. (|])-(|6|) for 
the solvent. This implies for the critical correlator the asymptotic law 
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CV\t) = // + ^(t/to)^ + O((t/t )- 2a ) ; a = 0. (16a) 

The nonergodicity parameter of the glass state, fj = C^\t — > oo), exhibits the y/a- 
singularity 



f j = f? + h j y/cx/(l-\) + 0(cT) , a^0+ , (16b) 
and the a-process initial decay is given by von Schweidler's law for t > t a and o — > 0— : 

c^\t) = ff\\- {t/fiY + o^t/fim , fi = (fi/h 3 y/ b c . (i6c) 

Let us introduce two ad hoc time scales for the description of the liquid relaxation outside 
the transient regime. The center of the /3-relaxation process, ri, shall be defined as the 
time, where the correlator has decayed to the plateau: C^'(ri) = The center of the a- 
process shall be defined as the time, where the correlator has decayed to 50% of the plateau: 
C^{ri) = f?/2. Some values are listed in Tab. |T|, and open squares and circles mark 
these a- and /3-relaxation times, respectively, in Fig. [|. The slowing down of the dynamics 
upon approaching the glass-transition point is reflected by the increase of the time scales 
with decreasing distance parameter |e|. The two step scenario emerges, because the ratio 
of the scales t^/tL increases as well. The a-decay leads to the a-peaks of the susceptibility 
spectrum, which are separated from the microscopic excitation peaks by a susceptibility 
minimum, as is demonstrated in Fig. ||. There, the a-peak-maximum positions, uj J max , and 
the minimum positions, w 3 min , decrease for — ► <p c — . The open squares and circles in Fig. 
d demonstrate, that co>4 ax ~ and uj J min ps 1/t^ as |e| — > 0. The two-step scenario implies 
that the ratio ^ax/^mm a ^ so decreases upon approaching the glass-transition point. Thus, 
for if — > <f c —, the a-peak gets more and more separated from the rest of the spectrum. In 
this limit, the plateau height is the relative area under the x versus- log u curve [|L6| : 

/logo;-' . poo 
mm x ij), '(u) dlogo;/ / dlogo;. (17) 

-oo J — oo 

Figure |3] demonstrates that for t ^ 3 the dynamics deals with oscillatory motion, i.e. 
with rotations and librations which are influenced by steric hindrance affects. If these effects 
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would lead to some fast decay towards the correlator's long-time limit, one would find a 
white- noise low-frequency fluctuation spectrum: C^"(oS) ~ C^"(u = 0). Equivalently, 
one would obtain a regular low-frequency susceptibility spectrum varying linearly with u, 
X^'^uj) oc lo, as is indicated schematically by the straight dashed-dotted line in the upper 
left panel of Fig. [5[ A linear susceptibility spectrum is obtained for the glass-spectra for 
uj <C l/t a , since the correlators approach the limit fj exponentially for t ^> t a . This is shown 
by the (e > 0) spectra in Fig. ||. Such regular spectra are also found for the low-frequency 
wings of the a-peaks, since the liquid correlators approach zero exponentially for t ^> r 3 a . 
At the bifurcation point, however, the critical decay leads to a power law spectrum which, 
according to Eq. (|16a|) , reads 

X (J) "(uj) = hj sin(a7r/2) T(l - a) (ut ) a + 0((c^ o ) (2a) ). (18) 

For t < tj, the correlators follow the critical decay decay if |<r| is small. Therefore the 
spectra are approaching the asymptotic dAlaw for l/t a <C uj <C l/t , as is demonstrated 
for the (x = 4) results in Fig. [|. The stretching of the first relaxation step leads to the 
strong enhancement of the intensity of the spectral minimum Xmm = X (^min) relative to 
any possible estimation of a white-noise-background spectrum. This enhancement also is 
exhibited by the experimental data reproduced in Fig. |I[ 

Let us consider the probability density P(r),t) = (S(r](t) — 77)) for the molecule's axis 
e(t) to have the projection rj(t) on its initial direction e: r](t) = e(t)e. Since S(rj(t) — rf) = 

1/2 + E"i (3 + 1/2) PMPAvit)), one gets 

00 

P( V ,t) = l/2 + J2 (j + 1/2) Pj(v) C^it). (19a) 

Thus, knowledge of the set of C^'(t), j = 1,2, is equivalent to knowing P(r),t). If the 
summation over j is understood with the cutoff l co , Eq. (|19a|) describes the evolution of the 
distribution with the initial value 

' c o 

P( V , t = 0) = l/2 + J2 U + 1/2) PM- (l 9b ) 
i=i 
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Figure ^ exhibits results for the small distance parameter — e = (ip c — tp)/ip c = 0.001 cor- 
responding to x — 3. The dotted lines exhibit P(i],t = 0)/10, calculated with l co = 7 for 
( = 0.80, and l co = 5 for £ = 0.33, respectively. Within the dynamical window, where the 
leading-order result for the /3-relaxation, Eq. ([15]), applies, one gets 

P( V ,t)=P c ( V ) + H( V )G(t) , (20a) 

oo oo 

P C ( V ) = 1/2 + £ V + V2) /; PM ; H(V) =E^' + X /2) fc, Pi(i7) • (20b) 

Thus the distribution relaxes towards the distribution P c (r]), which is frozen for tp = (p c . 
The relaxation does not exhibit any correlation between changes in time described by G(t), 
and variations with angle described by H{rf). This is the scenario expected for relaxation 
due to dephasing in the random distribution of sizes and shapes of the cages producing steric 
hindrance for the rotations. For ( = 0.80, the f3 regime extends from t = 10 to about 10 4 
as shown in Fig. |], and the upper panel of Fig. ^| exhibits the described phenomena for 
t = 10 2 and t = 10 4 . The /3-relaxation window is somewhat smaller for ( = 0.33, as will be 
discussed in quantitative detail below in connection with Fig. 0. The dephasing relaxation 
for this case is demonstrated in Fig. [J] for t = 10 2 and 10 3 . 

The beginning of the a-relaxation process follows von Schweidler's law, Eq. (|16cj) . It is 
identical with the end of the /^-process, and thus it is described within the scenario based 
on Eqs. (p0|). The most drastic difference between large- and small-elongation relaxation 
shows up for the a-process outside the von-Schweidler regime. For ( = 0.80, the probability 
decreases monotonically if the angle of the axis increases from its initial value O = 
to O = 7r. This is shown in the upper panel of Fig. ^ for t > 5 • 10 5 . As time increases, 
the probability for rj w 1 decreases, while it increases for 77 ~ — 1. Thus, the relaxation 
towards the equilibrium distribution P(rj, t = 00) = 1/2 is similar to what one would expect 
for diffusion on a sphere. For ( = 0.33, the correlators for odd j decay faster than the 
corresponding correlators with the even index (j + 1). This is demonstrated in Fig. |] and 
by the numbers in Tab. |T|. Therefore the a-process consists of an intermediate time 
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step leading to a probability distribution which is nearly symmetric with respect to the 
equator rj = 0. Only at later times, the symmetric distribution relaxes to the equilibrium 
one. Figure ^| shows that, already for the rather short time t = 10 2 , P(rj,t) exhibits a 
minimum. For t = 10 3 , there is an overshooting effect of the probability for e(t) = —e: 
P(j] = —l,t = 10 3 ) > 0.5; and this effect increases if the time increases to t — 10 5 . Thus, 
the relaxation pattern is that expected for a random process of large-angle flips of the 
molecule's axis. 



B. Dipole-versus-quadrupole relaxation for strong steric hindrance 

The equations for the nonergodicity parameters |10| imply that the fj increase towards 



unity if the coupling coefficients in Eq. (|T2|) are increased towards infinity. For this strong- 
coupling limit, one derives from Eq. (|19a|) that P(j],t) — > 5(r] — 1). Because of continuity, 
for strong steric hindrance and for t < t a , P(rj, t) is a narrowly-peaked distribution centered 
around rj ~ 1. Thus, one expects the expansion coefficients fj for not too large values of j 
to form a smoothly decreasing sequence of j: fi > f 2 > fj ~ (fj-i + fj+i)/2- Table | 
demonstrates this result quantitatively for ( = 0.80 and (p = ip c : 

/i C >/2 >/s >/!> lar geC (21a) 

In particular, the ratio (/f//!) of the relative strengths of the a-peaks for the dipole re- 
laxation, fi, and for the quadrupole relaxation, f%, is larger than unity. One cannot con- 
clude quantitatively from /f//! the ratio XiC^max) /x^C^max) °^ ^ e a_ P ea k heights, since 
the shapes of the spectra depend on j. However, Fig. ^ demonstrates that the two ra- 
tios are close to each other. One can also characterize the a-peak height relative to the 
microscopic-peak height, rj = X i u L.ax) I 'x^" '(^mic) ' or relative to the minimum intensity, 



3 A V max// A \ mic/ 

r 'j — X^" ( w max) / 'x^" '(^min) • F rom Fig- H) one infers r\jri ~ 4, and r[/r' 2 w 3, i.e. the 



(j=l)-versus-(j=2) enhancement effect appears even more pronounced. 

According to Eq. ( |16b|) , the nonergodicity parameters increase with increasing (ip — ip c ). 
On the other hand, 1 — /? > fj — /?. Therefore, hj must decrease if /? increases, so that 
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the strongly coupled parameters fj leave the asymptotic regime for Eq. ( |16b|) for similar 
magnitudes of a. Table | quantifies this result for ( = 0.80. In particular 

h 1 <h 2 , large (. (21b) 

The reasoning assumes fj to be large, and thus it cannot be applied for too large j. There 
is some jo, so that hj decreases with increasing j for j > jo- Within the frequency window, 
where the leading-order asymptotic law for the critical decay is valid, Eq. (fl8|) , one derives 
an enhancement of the (j=2) spectrum relative to the (j=l) spectrum, since: 

rtl/rt) = h 2 /h x . l/t a « u « 1/to . (22) 

The dotted lines c\ and c 2 in the upper panels of Fig. | demonstrate this result. 

For a strongly near-(?]=l)-peaked probability distribution P(t], t), one can approximately 
replace averages of functions of rj by the functions of the average (77). Thus, Lebon et al. 
concluded: fj = Pj(fi) p8"| . Specializing to <p — ip c , one quantifies the sequence of /? in 
terms of its first value : 

fj = Pj(fi) , C - 00 . (23a) 



Substituting into Eq. Q16b[) and specializing to a — > 0+, one can also quantify the sequence 



of hj by the first term hi'. 

hj = P'j{f°) hi, C - 00 . (23b) 

From Tab. | one infers, that for ( = 0.80 the error of Eq. ( |23a| ) for j = 2 (3,4) is as small 
as 0.1% (3%, 7%), and Eq. flZgEp is obeyed for j = 2 (3,4) within 5% (22%, 45%). 

The strong nonlinear couplings of the structural-relaxation modes require that all corre- 
lators enter the first relaxation step, the second relaxation step, and the equilibrium state 
nearly at the same respective time. This is demonstrated in Fig. ^ for ( = 0.80. The most 
striking manifestation of the coupling effect occurs at the center of the /3-relaxation window 
for er — > 0—. In this case, the factorization theorem, Eq. fli"5|), is valid. All correlators cross 
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their plateau at the same time, say Tp, where Tp is the zero of the j3 correlator G(t). Because 
of the scaling law, Eq. (5), one gets the result Tp = t_t a , i.e. 

r| = Lt a , a -> 0_ . (24) 

Here, t a is the scale from Eq. (|5a|), and i_ is the zero of the master function: g_(t-) = 0. 
For the HSS it reads i_ = 0.704 ||17|| . The open and full circles in Fig. § show that the 
asymptotic Eq. (|4]) is obeyed very well for ( = 0.80. Since the a-processes of C^\t) and 
C {2 \t) start at the same time t-t a and reach zero nearly at the same time, one expects from 
C^(Lt a ) = /{ > / 2 C = C^(Lt a ) that the decay time for C« is larger than that for 

r 1 a >rl ) large C (25) 

Furthermore, the (t) //f- versus- log t plot is somewhat steeper than the corresponding 
graph for j=2. This means that the stretching is larger for the (j=2) a process than 
for the (j=l) ol process. If one interpolates the decay functions by a Kohlrausch law, 
C {t)/ fj ~ ex P [ — {t/ T a)^ j ] 5 the stretching exponent for j=l is larger than that for j=2: 

(3i>(3 2 , large (. (26) 

Stretching can also be quantified by the width w at half height of the a-peak of the suscep- 
tibility spectrum. For ( = 0.80 our model yields for j = 1(2,3,4) w = 1.16(1.25, 1.37, 1.50) 
decades. The Kohlrausch processes leading to the same w require stretching exponents 
P = 0.99(0.90,0.82,0.74). 

The derivation of the inequality for the time scales can be put on a quantitative level by 
combining Eq. (|16cj ) with the two inequalities in Eqs. (pi]). One gets in analogy to Eq. (p5|): 
fl > f^. The a-relaxation law for the C {i) (t) holds in analogy to Eq. (0): C {j \t) = Cj(t). 
If the shape function Cj{i)/ would be independent of j, the ratio t*/t^ would be equal to 
the ratio f^/fj = [/i/i2/^i/|] 1 ^- But the latter is about 2.1 times larger than r^/r^. 
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C. Dipole-versus-quadrupole relaxation for weak steric hindrance 



There are two universal phenomena which are relevant for a discussion of the dynamics for 
weak steric hindrance. The first one concerns the limit ( = of the center-of-mass correlator 
$(j = ji = q, t) = & s q (t), which is identical to the tagged-particle-density correlator of 
the larger of the two spheres forming the dumbbell. If the radius of this sphere, say di, is of 
the same order or larger than the radius d of the solvent spheres, the steric hindrance is very 
effective. In this case, $f(i) exhibits the canonical bifurcation scenario if (p crosses (p c , as was 
discussed comprehensively in Ref. [IB]. This implies that for di > d the (j = 0) correlators 
exhibit only a smooth ^-dependence for £ decreasing to zero. A side remark shall be added 
to this conclusion. If the ratio of the diameters di/d of a sphere moving in a glass of hard 
spheres decreases towards zero, there occurs a percolation transition at some critical value 
(d\/d). This is a type-A transition, i.e. a bifurcation where the Lamb-Mossbauer factor 
decreases continuously to zero for (d\ — df) approaching zero from above |I5|JI^] . Because 
of continuity, it is obvious that for a dumbbell built of sufficiently small spheres, d\ < d\, 
there will be a type-A transition if the elongation ( decreases to some critical value (* > 0. 
If ( crosses (*, the dynamics changes from one dealing with molecules localized in the hard- 
sphere glass to one dealing with delocalized molecular motion. This small-C phenomenon 
for small molecules is not considered in this paper. 

The second universal phenomenon deals with a type-A transition resulting from the 
fact, that for top-down symmetrical molecules the MCT equations of motion of the even- 
j correlators decouple from the odd-j ones |To|j26|] . The even-j correlators couple to the 
function Q s q {t), and thus the conventional transition scenario of this correlator enforces the 
same for all other correlators with even j. However, such coupling does not exist for odd 
j. For large (, this results in no considerable effect. But all coupling coefficients in the 
equations of motion approach zero for odd j if £ tends to zero. Consequently, for all if > 
ip c there is some critical elongation ( c ((p) for a type-A transition. For the studied model 



C c (<£>) < Cc(y?c) = 0.296 [T0JT1]. Choosing ( sufficiently close to (d^), it can happen that 
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for odd j: fj < fj + \ or even fj < fj + 3 [10|. The transition at Cc(y?) shall not be studied in 
this paper. For the demonstration of the small-(C — Cc(Vc)) phenomena, the value ( = 0.33 
has been chosen so large that the canonical sequence for the plateau values in the (q = 0) 



limit, Eq. ( |21a[ ), is not violated, as is quantified in Tab. |I[ But it is chosen so small, that 



the precursor effects of the type-A transition seriously influence the results for the dumbbell 
dynamics. Thereby, the results are also representative for such cases, where the type-A 



transition singularity is avoided [15] due to a weak breaking of the top-down symmetry of 
the solvent-solute interaction. 

The even-j correlators show the conventional behavior. Therefore, the discussion of their 
trends with decreasing ( for fixed j can be held analogously to that given in Sec. IIIB for the 
trends with increasing j for fixed large (. Thus one understands that the (j=2) a-process for 
( = 0.33 is weaker, faster, and more stretched than that for ( = 0.80, as it is demonstrated 
in Figs. f| and |5|. The halfwidth of the a-peak for j = 2(4) is w = 1.66(1.86) decades 
as for a Kohlrausch process with exponent (3 = 0.67(0.59). Notice in particular from Fig. 
[|, that the /3-relaxation scale r| for the x = 3-result is close to the (- and j-independent 
number t_t CT from Eq. (]2~4]). For x = 2, the asymptotic formula is obeyed reasonably, but 



the preasymptotic corrections are larger for ( = 0.33 than for ( = 0.80. 

The most obvious precursor of the type-A transition is the suppression of the plateau 
values fj for odd j. This leads to a violation of the rule (/i + fa)/ 2 ~ fa, as is quantified in 
Tab. |. The general qualitative reasoning from Sec. IIIB explains, that the suppression of 
fi is connected with an enhancement of h\. hi(( = 0.33)//ti(£ = 0.80) m 15. The amplitude 
h% is given by the resolvent of the so-called stability matrix, and at a type-A transition the 
resolvent exhibits a pole [|l7||16|l . Hence hi(( — > Cc(y))/^2(C ~" ¥ Cc(y?)) - > oo, and the regular 



trend, Eq. ([21 b|) , is reversed: 

hi > h 2 , small C- (27) 

For our example one infers from Tab. | that hx/hi ~ 4.2. According to Eq. (P2"D, the 
critical spectrum for the dipole relaxation is considerably larger than that for the quadrupole 
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relaxation, as is demonstrated by the dotted lines in the lower two panels of Fig. [|. 

Combining Eq. fl2~T|) with von Schweidler's law, Eq. ( |16c| ), one concludes that the C^(t)- 
versus-log t curve crosses its plateau f{ much steeper than the C^- 2 ' (t)-versus-log t curve. This 
is illustrated in the lower panel of Fig. |j. Hence the a-relaxation of the (j=l) response is 
faster than the one of the (j=2) response: 

rl<rl, small (. (28) 

Again the order for large (, Eq. (p5|), is reversed. From Tab. ||, one infers for x = 3: 
T a/ T a = 0.12. Accordingly, the a-peak positions for the (x = 3) spectra for j—1 and j=2 
in the lower panels of Fig. [5] differ by about one order of magnitude. For the ratio of 
the von Schweidler scales in Eq. ( |16cj ), one gets f^/f 2 = [h 2 fi /hif^] 1 ^ — > for £ — ► ( c , 
and this identifies the smallness of the ratio r 1 /r 2 as a precursor of the type-A transition. 
The preceding discussion is valid more generally and explains that all the odd-j correlators 
decay faster than the nearby even-j ones. As a result, the probability distribution P(j],t) 
approaches first an even function of r], before the equilibrium value 0.5 is obtained, as is 
demonstrated in the lower panel of Fig. |[ 



D. /3-relaxation scaling 



The factorization theorem for the /3-relaxation, Eq. (jig) , means that the rescaled corre- 
lators c^'(t) = (C^'(t) — fj)/hj are given independently from j by the /3-correlator G(t) of 
the solvent. The latter obeys the scaling law, specified by Eqs. (^)-(|6a|). For fixed rescaled 



time t = t/tfj, the cited formulas deal with the results correctly up to order y \a\ [0. The 



leading corrections are of order |er|, and they explain the range of validity of the leading 



results for separations e |T7[| . Figure demonstrates these statements. On a 10%-accuracy 
level the leading-order results describe 4%(18%; 45%; 20%) of the decay of the correlators 
around the plateau for ( = 0.80, j = 1(C = 0.80, j = 2; C = 0.33, j = 1; C = 0.33, j = 2). 
These decay intervals are indicated in Fig. [| by vertical lines. For e = —0.001, the corre- 
sponding dynamical window extends from about t = 10 to about 10 5 , while it extends from 
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about t = 3 to about t = 100 for e = —0.01. This discussion requires a reservation: The 
corrections to the scaling results can lead to such a violation of Eq. (|24|), which appears as 
an offset of the plateau |T7|]. This offset can be noted in the lower panel of Fig. [7] for the 
odd-j results. The good description of the /3-decay of the (( = 0.33) results for odd j holds 
only after a correction of the offset. For t > 10 4 , the correction effects cause the c {j \t) for 
C = 0.80 to differ from G(t)\ one infers from Fig. [7] that the c^\t) increase with increasing j. 
The general results for the theory of the corrections imply, that then also the c^' (t) increase 
with j for t < 10 jl7| . The (t)-versus-log t curves do not intersect for i-t a but they touch 



each other as is demonstrated in the upper panel of Fig. [7]. Corresponding results also hold 
for ( = 0.33 after the mentioned offset is eliminated. 

Equations (H) and (|T5|) lead to the scaling law for the susceptibility spectra: x i u ) I hj = 



c o-X±(u^o-)- The master spectra x±(u>) = &g±(&) are given by the Fourier-cosine transform 
g±(uj) of the master functions g±(i). The master spectrum for the glass state describes the 
crossover from a regular spectrum for small rescaled frequencies, x+(ci> ^ 1) ^ to the 
critical spectrum at large rescaled frequencies, x+{^° ^ 1) °^ & a - It deals with the knee 
exhibited by the spectra for e > and x = 3, 4 in Fig. [5]. The master spectrum for the 
liquid describes the crossover from the von-Schweidler high-frequency tail of the a-peak, 
X-(lj <^ 1) oc l/uj b , to the critical decay for large rescaled frequencies, X-(& 3> 1) oc Cj a . 
The results describe in the small-a limit the /3-relaxation minimum as it can be seen in 
Fig. H for the x = 3 and x = 4 results. In particular, the factorization theorem explains 
why the spectral minima c<Vin are located at the same position independently of j and £• 
The leading-order formulas imply co> min = cj min /t CT , where cu min denotes the minimum of the 
master spectrum x~- F° r the hard-sphere system, one gets u min = 1.56 [ |I~7| . 



Obviously, the /^-relaxation scaling laws can describe the susceptibility minimum only for 
such small distance parameters, for which kVi n is located in that frequency window where 
the (a = 0) spectrum exhibits the asymptotic u; a -law, Eq. (|T8|). Figure [| shows that for the 
model under study this window is restricted to uo < 0.01. This means, that u; miri has to be 
located about three decades below the peak of the microscopic susceptibility spectrum. For 
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to > 0.01, the critical spectrum is modified by crossover effects to the transient dynamics. 
The susceptibility minimum with > 0.01 is due to the crossover of the a-peak tail to the 
microscopic excitation spectrum; it cannot be discussed by the universal asymptotic laws 
for the MCT bifurcation. One concludes from Fig. [|, that |e| < 10~ 2 needs to be satisfied 
in order to apply the /3-scaling laws for the model under study. 



E. o-relaxation scaling 

The a-relaxation scaling law reads for the reorient at ional correlators in analogy to Eq. 

(0): 

CW(t) = Q(f) I t = t/l. (29) 

The e-independent master function Cj obeys as initial condition the von Schweidler law: 
CU'yt) = fj — hjP + 0(t 2b ). The superposition principle for the susceptibility spectra 
reads x^"^) = X {&) with uj = cut' a denoting the rescaled frequency. The e-independent 
master spectrum is given by the Fourier-cosine transform of the master correlators x (ui) = 
uC^"(u). Consequently, the above defined a- relaxation time scales r 3 a and susceptibility- 
maximum positions u; max read 

4 = H > <ax = , (so) 

where P is defined by C^'(P) = f?/2 and Cji denotes the peak frequency of X U) "(Co). The 
scaling law implies that a representation of C^\t) as a function of the rescaled time t/r 3 a 
should superimpose correlators for different distance parameters e on the common curve 
C ij) (t/t j ). Asymptotic validity means that the log(£/r^) interval, where the scaling law is 
obeyed, expands to arbitrary size for e — >• 0. A corresponding statement holds for the repre- 
sentation of the susceptibility peaks as functions of the rescaled frequency. The corrections 
to the leading-order asymptotic laws are the larger, the larger the critical amplitude hj is 
[ r^JIBfl . Figures | and |3| demonstrate, that the described scenario for the evolution of the 
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a-process is valid for ( = 0.80, and also for ( = 0.33 in the case j=2. For strong steric 
hindrance, the a-scaling law works for larger values of ((p c — ip), than the /3-scaling law. This 
is so, because the leading corrections to the a-scaling law are of relative size (9(|e|), while 
they are of relative size (9(-\/je[) for the /9-scaling law |T7| . 

Figure § demonstrates a drastic (j=l)-versus-(j=2) effect of the a-scaling for ( = 0.33. 
The dipole correlators do not exhibit the superposition principle for |e| > 10 -4 , nor do the 
correlators for the other odd values of j. For j = 1 the plateau emerges only for extremely 
small values of the distance parameter |e| < 10~ 4 . The a-peak heights of the dipole spectra 
decrease with decreasing |e| in Fig. [5] in contradiction to the scaling-law prediction. This 
anomaly is caused by the large size of the critical amplitude h\, which was explained in 
connection with Eq. ([271). More precisely, it is caused by the large percentage of the decay 
of C^(t) described by the /3-scaling law as is indicated by the vertical lines in Fig. |j. To 
formulate this observation quantitatively, let us remember that the decay of the correlator 
near the plateau is described by Eqs. (Q), and (p~5|) : C^(t) = /? + hj^J\a\g_(t/t a ). The 
master function g_(t) for small positive values and all negative ones is well approximated 
by: g{t) = —Bt b + Bi/(Bt b ). Here, B\ is determined by the exponent parameter A and for 
our solvent model reads B\ = 0.431 |T7J. Thus one gets for C^\t) < fj within the window 



for the validity of the /3-relaxation scaling law: 

C (j) (t) = /? - v ^\h J {B(t/t (7 ) b - B l /[B(t/t a ) b }} . (31a) 

The leading corrections to this formula can explain the possible offset of or, equivalently, 
of the scales t a [|17]], which was noticed above in connection with Fig. |7| for ( = 0.33. 
Equation ( |31a| ) can be rewritten as C^\t) = [/? — hjV 1 ] + hj\a\Bi/i b . Here, the bracket is 



the a-scaling-law description of the initial part of the a-process, and the term proportional 
to Bi is the leading correction. The correction term to the a-scaling law deals with that 
part of the /3-process below the plateau, which is not given by the von Schweidler's large-t 
asymptote. Therefore, one can write the for the a-process for not too large values of rescaled 
time t: 

24 



C®(t) = C {j \t) + hjlalBj? . (31b) 
The analogue of this formula was shown in Ref. [17] to describe the evolution of the a-process 



of the density correlators of the HSS perfectly for |e| < 0.1. It was also shown, that the 
corresponding spectrum describes the susceptibility peak to increase above the scaling-law 
constant ^"(^max) if the separation |e| increases from 10~ 2 to 10 _1 . 

In Fig. ||, the evolution of the (j=l) a-process for small steric hindrance is reexamined. 
Instead of rescaling the time with the ad-hoc scale r^, the theoretically motivated scale t' a is 
chosen. One recognizes, that the found scenario does not exhibit any qualitative peculiarity 
anymore, compared to what is presented in Fig. [8] for Q = 0.33 and j=2. The (j=l)-versus- 
(j=2) anomaly is identified as an anomaly of the size of the corrections only. In the case 
of the small elongation, the distance parameter |e| has to be taken almost two orders of 
magnitude smaller for j=l in order to render the corrections to the a scaling as small as 
found for j=2. For |e| ^ 10~ 3 , even including the leading corrections to von Schweidler's 
law, one can explain the relaxation from the plateau only up to some offset in the time scale. 
This is demonstrated in Fig. |9] by the dotted lines for x = 2,3. 

Two remarks concerning tests of the second scaling law shall be added. The definition 
of the time scale used in Fig. |] was arbitrary. Let us consider more general definitions 
to be parameterized by a number k > 1 and denoted as r k . The subscripts a and j shall 
be dropped for the sake of simplicity, and the definition shall be: C^(rfc) = ff/k. If the 
scaling law is valid, one finds in analogy to Eq. (p0|): r k = t k t' a . Here, i k is defined by 
C®{i k ) = fpk. In this case, the choice of k is irrelevant, since the ratio of two different 
scales is e-independent: r kl /r k2 = t^/t^- However, if preasymptotic corrections are present, 
the scales are not equivalent. The range of validity of the superposition principle expands 
from large to small rescaled times. This follows from Eq. (|31b| ) and is demonstrated in Fig. 

One gets for k\ < k 2 : 

Tki/hi < T k 2 /h 2 < t'a ( 32 ) 

For a detection of the superposition principle for an as large as possible value of |e|, one should 
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therefore choose an as large as possible value of k for the rescaling procedure. Thereby, the 
artificial crossing point of the rescaled curves at t — is suppressed as much as possible. 
Otherwise, one introduces a time scale for the characterization of a decay process which 
cannot be characterized by a single scale. The outcome of this ill-defined procedure is 
demonstrated in the lower left panel in Fig. [S[ In this case, r\ = r^ =2 is a parameter 
extracted from the correlator which, according to Fig. [7|, is adequately to be specified by 
the two scales c a and t a of the /3-relaxation scaling law. The dashed line in the inset of 
Fig. ^| demonstrates explicitly, that the scale does not exhibit the asymptotic behavior 
for |e| > 10~ 3 . The asymptotic law t\ = t-^! a is followed only for e < 10~ 4 . 

The second remark concerns the determination of the exponent 7 entering the power-law 
behavior for the a-relaxation time scale, as specified by Eqs. (|6l5| ) and fl30l). These results 
are based on the validity of the scaling law [[UJ. Therefore, one cannot appeal to MCT if 
one fits power laws for scaling times for cases where the scaling law is violated. The dashed 
line in the inset in Fig. |9| demonstrates, that the scale for |e| > 10~ 3 can be fitted well 
by a power law for a two-decade variation of the distance parameter |e|. The identified 
effective exponent 7' < 7 describes the variation of t\ over three orders of magnitude; but 
nevertheless 7' has no well defined meaning for the discussion of our model. 

IV. CONCLUSIONS 

Solving the MCT equations of motion for the dynamics of a hard-sphere dumbbell moving 
in a hard-sphere liquid, first-principle results have been obtained for the evolution of the 
glassy dynamics of the reorient at ional degrees of freedom of a molecule. It was found that 
one has to distinguish between two scenarios, namely between strong steric hindrance as 
found for large elongations ( of the dumbbell, and weak steric hindrance as found for small 
elongations. 

For strong steric hindrance, the mode-coupling coefficients for the reorient at ional degrees 
of freedom in Eq. flT2]) are of the same order as the ones entering Eq. (|2]) for the description 
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of the translational degrees of freedom of the solvent. The dependence of the various pa- 
rameters on the angular-momentum index j is similar to the dependence on the wave vector 
q. One has to view j(j + 1) as the analogue of q 2 . While the g-dependence reflects the de- 
composition of the direct solute-solvent correlations in plane waves the j-dependence reflects 
the decomposition in spherical harmonics. Hence one finds, that all reported results on the 
j-dependence of the reorientational correlators C^'(t) are similar — and can be explained 
in a similar manner — as known from the previous work on the tagged-particle-density 
correlators in simple liquids (T^j2^|. In particular, it was shown that with increasing 

j the a-peak-strength parameters /?, Eq. ( |21a| ), the a-relaxation-time scales r^, Eq. (^5j), 
and the stretching exponents (3j, Eq. fl26|), decrease. These findings reproduce qualitatively 
the three general differences between dielectric-loss and depolarized-light-scattering spectra, 
which were discussed in Sec. I in connection with Fig. [I]. Because of Eq. (|la|), the re- 
laxation of the correlator $ follows that of the kernel m. Therefore, the a-relaxation time 
scale of the (q — * 0) density fluctuations, say r°, is larger than the corresponding scale of 
the longitudinal elastic modulus m q=0 (z), say t™. For strong steric hindrance, the decay 
of the cage is the prerequisite for the reorientation of the molecule, and therefore r° < t\. 
Thus one expects the fourth general feature of the a relaxation listed in the introduction: 
r^/r™ > 1. For our model one gets for e = -0.01: t™ = 130, r° = 240, t 2 = 920. The ratio 
T a/ T a ~ 7 is of the same order as cited in Sec. I for PC and Salol. 

The mode-coupling coefficients in Eq. ([12]) decrease to zero if j tends to infinity. Thus the 
solutions for large j are sums of many small terms, which are not strongly correlated. Each 
term exhibits the short-time von Schweidler-law behavior for the a-relaxation: C^'(t) — /? oc 
(t/t' a ) b . Therefore, one expects for C^'(t) the characteristic function of the stable Levy 
distribution, exp [(— tFj) b ] [[|. For the density correlators of the solvent, Fuchs has worked 
out the limit behavior for the a-relaxation master function for q — > oo and showed how the 



Kohlrausch law with (3 — b arises |3(J. We suspect that a similar derivation can be done for 
the reorientational correlators. Therefore, we conjecture that the sequence of Kohlrausch 
exponents (3\ > fa > fa... converges towards the von Schweidler exponent b. Molecular- 
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dynamics-simulation data for a model of water have been interpreted consistently within 
the standard MCT scenario pT|-|33|j. In particular, the C^\t) exhibit the conventional 



behavior p% . One concludes that water exhibits strong-steric-hindrance effects. Therefore 
it is reinsuring that the sequence of the first five Kohlrausch exponents (3j decreases with 
increasing j monotonously towards the von Schweidler exponent [|35|| . A further general 
result, namely the increase in the initial part of the series of critical amplitudes, Eq. (|21b|) , 
is also found in the simulation data for j < 3 |35] . 



Figure p] exhibits as full lines the (j=l) and (j=2) spectra calculated for ( = 0.80. The 
lines for x = 1 and x = 2 are the ones discussed in Fig. ||, and the other two refer to x = 1.33 
and x = 1.67, respectively. In order to transfer the MCT results, which are calculated with 
ad hoc units specified in Sec. II, to the units used by the experimentalists, one has to 
introduce three scales. The first and second scale transfer the calculated dimensionless 
normalized spectra x i u ) f° r J = l an d J=2 to the units used by the experimentalists for 
their dielectric-loss and depolarized-light-scattering spectra, respectively. The third scale 
shifts our frequency scale to the GHz scale. In the double-logarithmic representation, the 
first two scales define an overall vertical shift of the diagrams in Fig. [|, while the third 
scale defines a horizontal shift of the figures. Intending to compare data for PC for different 
temperatures with the MCT results for different packing fractions if, one gets a mapping of 
the T-scale onto the <f -scale via Eq. @ . The result is shown as an inset in Fig. p]. The inset 
also includes the point with coordinates of the critical packing fraction of the solvent ip c and 
the critical temperature T c = 180K. This value for T c was determined for PC by analyzing 



neutron scattering data |36|] , and has recently been corroborated in an MCT analysis of 
various other PC experiments |57j]. Our results in Fig. [I] describe the evolution of the 
two types of PC spectra semi-quantitatively. In particular, the extrapolation of the T-ip- 
relation yields a reasonable estimation of the critical temperature for that system, which 
is demonstrated through the dashed line in the inset. There is no obvious reason why the 
studied dilute solution of hard symmetric dumbbells in a hard-sphere solvent should produce 
spectra, which are similar to the data for PC. We consider the found similarities to a large 
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extend as accidental. The theoretical curves are added in Fig. |] with the mere intention 
to justify the conclusion: the model studied in this paper and our choice of parameters are 
relevant for achieving an understanding of experiments in glass-forming molecular liquids. 

In order to further corroborate the preceding conclusion, let us consider Fig. [10]. The data 
points exhibit a susceptibility spectrum of PC measured by incoherent-neutron-scattering 
spectroscopy for the wave vector q = 1.3A 1 j38|. A remarkable feature of the a-peak 
spectrum compared to the spectra shown in Fig. [I] is that it is less pronounced relative to the 
spectrum of the microscopic excitation band and that is is more stretched. The two dashed 



lines in Fig. |10| exhibit the spectra for the center-of-mass correlator = &(q,j = 0,/i = 
0, t) for q = 7.4 in order to emphasize, that this leading approximation for the scattering 
function cannot easily explain the experimental findings. The shown a-peaks of Q s have a 
half width of w — 1.34 decades as produced by a Kohlrausch process with exponent f3 = 0.84. 
The scattering function F q (t) is a sum over the contributions of the molecule's constituents 
and hence it is a superposition of the density correlators for all angular momentum indices 
j. For the symmetric dumbbell one gets up to some normalization constant [TO 



F q (t) = (2j + 1) bj(q C/2) 2 $(q,j,0,t) , (33) 

3 

where bj(z) are the spherical Bessel functions. The full lines in Fig. 10 exhibit the spectra 
for F q (t) for q = 7.4. The a-peaks have a halfwidth of w — 1.61 decades as produced 
by a Kohlrausch law with stretching exponent (3 = 0.69. The frequency was rescaled as 
explained in connection with Fig. [I] and the scale for the spectra was adjusted to meet the 
one of the data. Comparison of the full line with the dashed one for x = 2 shows the features 
distinguishing the a-processes of F''{uj) from that of The result calculated for x = 1 

shows that the finding for our model semiquantitatively accounts for the a-peak data. 



Some side remarks considering the comparisons in Figs, p] and 10 might be useful. A 



schematic-model analysis of the PC data gave the exponent parameter A ~ 0.75 |37], in 



good agreement with the values found from analyses of the susceptibility minima with the 



/3-relaxation scaling laws |19|j20|j38|| . The value is close to the result A ~ 0.74 for the hard- 
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sphere system, Eq. This accident ensures that the master function for the susceptibility 
minimum and the values of all anomalous exponents of PC agree within the experimental 
uncertainties with the corresponding quantities of the model studied in this paper; and this 
is a prerequisite of a successful fit. Accidently, the ratio of the a- relaxation times r„/r^ 
noted in Tab. [TT] for ( = 0.8 is only a bit larger than the ratio of the a-peak-maximum 
positions of PC, ^max/^max) exhibited in Fig. 1; and this is another request for a reasonable 
fit. Since the ratio decreases with decreasing (, some ( < 0.80 could be chosen to reproduce 
the specified (j=l)-versus-(j=2) effect quantitatively. Nevertheless, it is remarkable that 
the fit in Fig. [l| reproduces the ratios of a-peak maximum intensity to /3-minimum intensity 
X {umaxj/x (k-Wn) reasonably well for both values of j. Neither it is trivial, that the 
model reproduces reasonably the (j=l)-versus-(j=2) effect for the stretching. 

A new liquid-glass-transition scenario is predicted which is referred to as the regime of 
weak steric hindrance for reorientational motion. It is characterized by (j=l)-versus-(j=2) 
effects, more generally by odd-j-versus-even-j effects, which are qualitatively different from 
the results described above as strong-steric-hindrance results. The new scenario occurs, if 
precursor phenomena of a type-A-transition between two nonergodic states strongly influ- 
ence the asymptotic results for the conventional MCT bifurcation. The scenario appears if 
the particle interactions deviate not too strongly from spherical symmetry, e.g., if a linear 
molecule exhibits only small deviations from a top-down symmetry and if there are not too 
large elongations. Six features characterize the weak-steric-hindrance scenario. First (i), the 
plateaus /J for the reorientational correlators for odd j are suppressed in comparison to what 
one would expect by interpolating or extrapolating the values for nearby even-j plateaus 
(Tab. |). Most importantly (ii), the critical amplitude hi is larger than h,2, Eq. (|27|), so 
that the canonical ordering of the hj for small j, Eq. ( |21b| ), is reversed. Third (iii), the 
percentage of the decay of the reorientational correlators C^\t), which can be explained by 
the leading- plus next-to-leading-order asymptotic formulas for the /3-relaxation is larger for 
j=l than for j=2; as is indicated by the vertical lines in Fig. 4. The structural relaxation of 
the reorientations is dominated by large-angle flips (iv), as shown in Fig. |^ for the dumbbell 

30 



with elongation ( = 0.33. The a-relaxation-time scale for j=l is smaller than for j=2, Eq. 

(v) so that the canonical order of the a-relaxation scales, Eq. (^), is reversed. This 
can cause the C^> (t)-versus-log t graphs for j—1 to cross the graphs for j—2, as is shown 
for the ( = 0.33 results in Figs. |3| and [|. Finally (vi), for distance parameters |e| > 10~ 3 , 
where the conventional C^(t)-versus-logt/(r^) plot exhibits the a-relaxation scaling law 
for j=2, the correlators for j—1 do not show the validity of the superposition principle, as 
is demonstrated in Fig. ||. Nor does the scale r*, defined as the time for a 50% decay of 
the a-relaxation correlator, exhibit the power law behavior with the correct exponent 7 as 
is shown in the inset of Fig. |j. 

A side remark concerning a molecular-dynamics study of the evolution of glassy dynamics 



in a Lennard- Jones-dumbbell liquid by Kammerer et al. p9|-[41]1 might be in order. It was 



reported that the correlators dealing with translational degrees of freedom and also for 
the ones for the reorientational dynamics for angular index j / 1 could be interpreted 
qualitatively within the universal asymptotic MCT formulas. However, the evolution of the 
dipole correlators did not fit into the standard MCT pattern. It was found that hi > h 2 and 
T a < T a- A drastic violation of the a-scaling law was noted quite similar to what is exhibited 
in the lower left panel of Fig. [5]. The scale t\ exhibited a deviation from the asymptotic law 
e~ 7 , but a fit by r a oc e -7 ' with 7' as discussed in the inset of Fig. || was possible. These 
simulation results for j=l differ from those for water simulations [j32H35"j as well as from the 



experimental findings for propylene carbonate quoted in Fig. [I]. However, they agree with 
the features (ii), (v), and (vi) specified in the preceding paragraph. Moreover, the property 
(iv) concerning the large angle flips is also obtained in Ref. |39| . Accidently, the x = 2 
results in the lower panel of Fig. [| show that the minimum position of the j=l spectrum 
exceeds that of the j=2 spectrum by nearly one order of magnitude, in agreement with the 
corresponding finding in Ref. ||41|| . Furthermore, the a-peak variation with x shown in the 
lower left panel of Fig. |5] is in qualitative agreement with that reported in Ref. ||39|| . In view 
of these observations it does not seem impossible, that the scenario studied in Refs. p^-fTI 



fits into the framework of the ideal MCT. However, it is not clear, whether or not the results 
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of Ref. |39|-f4T| can be explained by our theory for type-A precursors of a dilute solution of 
molecules. First, the simulation results for the dumbbell liquid do not exhibit a particular 
decrease of /f relative to /|. Second, the /3-relaxation scaling has not been documented for 
the dumbbell liquid and so it is unclear, whether or not the feature (iii) holds for that case. 

Summarizing, it shall be emphasized that all qualitative features for the evolution of the 
structural relaxation studied in this paper have been explained by means of the formulas for 
the leading-asymptotic expansions and their leading-order-correction formulas for the bifur- 
cation scenario. In this sense, these asymptotic formulas can be considered as the essence 
of MCT. However, in order to explain the characteristic (j=l)-versus-(j=2) differences for 
the relaxation patterns, it is necessary to also understand the general trends of the nonuni- 
versal parameters with variations of wave- vector q and angular- momentum index j. And 
this requires to use MCT as a microscopic theory based on the knowledge of the equilibrium 
structure. 
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TABLES 





TABLE I. Plateau values /? and critical 


amplitudes hj 








C = 0.80 




C = 0.33 


J 

1 

2 
3 

4 


h 
0.943 

0.835 

0.701 

0.540 


hj 

0.13 
0.35 
0.55 
0.68 


h 
0.303 

0.286 

0.052 

0.006 


h ■ 

1.94 
0.46 
0.46 
0.13 


TABLE II. Time scales and ri 






C = 0.80 


C 


= 0.33 


rl 
rl 
rl 
rl 


x = 2 
2.75 x 10 3 
9.21 x 10 2 

4.40 x 10 2 

2.41 x 10 2 


x = 3 
7.65 x 10 5 
2.56 x 10 5 
1.20 x 10 5 
6.43 x 10 4 


x = 2 
9.29 x 10 1 
1.85 x 10 2 
5.98 x 10 1 
6.55 x 10 1 


x = 3 
6.51 x 10 3 
5.39 x 10 4 
4.43 x 10 3 
1.74 x 10 4 


4 


1.37 x 10 1 

1.31 x 10 1 

1.32 x 10 1 
1.25 x 10 1 


5.42 x 10 2 
5.34 x 10 2 
5.32 x 10 2 
5.22 x 10 2 


3.33 x 10 1 
1.03 x 10 1 
2.53 x 10 1 
1.00 x 10 1 


8.11 x 10 2 
4.80 x 10 2 
7.24 x 10 2 
4.74 x 10 2 
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FIGURES 

FIG. 1. Susceptibility spectra x" °f propylene carbonate (PC, symbols) and solutions ob- 
tained for a symmetric hard-sphere dumbbell with elongation £ = 0.80 immersed in a hard-sphere 
solvent (full lines, see text for details). The symbols represent dielectric- loss spectra measured 



by Schneider et al. |Tj| (upper panel) and depolarized-light-scattering spectra of Du et al. [2C] 
(lower panel) for temperatures as indicated. The full lines are calculated for the distance param- 
eter e = {ip — ip c )/ip c = — W~ x with x = 1,1.33,1.67 and 2 for angular momentum index j=l 
and 2, respectively. Computed frequencies have been rescaled by a factor of 10 to meet the ex- 
perimental GHz scale. The calculated susceptibilities have been divided by 2.8 for the j=l case 
in order to normalize the spectra at lo/2tt = 2GHz. The inset exhibits packing fraction ip versus 
temperature T for which the spectra are fitted. Here, the critical value of the hard-sphere system, 
ip c = 0.516, corresponding to the critical temperature of PC, T c « 180K, was added. The dashed 
line demonstrates the extrapolation from the found <^-T-mapping to T c . 



FIG. 2. Angular dependent solute-solvent pair-distribution function g(f, fi), calculated within 
the Percus-Yevick theory, for a top-down symmetric solute molecule consisting of two equal fused 
hard spheres with elongation ( = 0.8 (upper panel) and C = 0.33 (lower panel). The shown x— z 
plane contains the molecule axis. Grey corresponds to g(f,Cl) ~ 1, dark and white areas show 
regions with higher and lower probability to find a solvent particle, respectively. The cut through 
the dumbbell is shown hatched. The diameter d of each sphere is chosen to match that of the 
surrounding solvent particles. The unit of length is chosen here and in all following figures such 
that d = 1. The packing fraction of the hard-sphere solvent is at the critical value ip c = 0.516. 
g(f, O) was approximated using a Legendre-polynomial expansion with angular momentum indices 
up to j = 16. 
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FIG. 3. Correlators $ for the wave vectors q = 7.0 and 10.6, elongations £ = 0.80 (upper two 
panels) and C = 0.33 (lower two panels), angular indices j = 0,1,2, and helicity index \i = as 
functions of the logarithm of the time t. The unit of time is chosen here and in all following figures 
such that the thermal velocity v of the solvent is unity. Correlators are shown as full lines for 
j = 0, 2 and as dashed lines for j=l. The solutions at the critical packing fraction are marked by a 
c and are shown in dotted. The packing fractions are parameterized as (tp — ip c )/ip c = e = ±10~ x , 
and x = 1,2,3,4 was chosen. Solutions for the glass states, e > 0, are only shown for £ = 0.80, 
j=l. Correlators are truncated where necessary to avoid overcrowding of the figure. 

FIG. 4. Reorientational correlators C 3 (t) for j=l and j=2 for the two elongations £ = 0.80 
and C = 0-33 as function of log 10 1. The solutions at the critical point are shown in dotted and are 
marked by Cj. The plateau values /? are marked by horizontal lines. The distance parameter is 
chosen as e = (ip — ip c )/ip c = — 10 _a: with x = 3 (slower decay) and x = 2 (faster decay). Open 
circles and open squares mark the characteristic time scales and for the a- and /3-process, 
respectively. The full circles and squares mark the time scales 0.704 t a , with t a from Eq. (|6a|) , 
and t' a from Eq. (|6~b|), respectively. The vertical lines indicate the decay interval described by the 
asymptotic formulas for the /3-process (see text, cf. Fig. ^). 

FIG. 5. Double-logarithmic presentation of the susceptibility spectra x^"^) = wC^'"(w) 
for angular-momentum indices j=l and j=2 for elongations £ = 0.80 and £ = 0.33. Spectra for 
the critical packing fraction ip = ip c are shown in dotted and are marked by Cj. The distance 
parameters are e = ±10~ x with x as given in the panels. In the upper left panel, a regular 
susceptibility spectrum, \" oc u>, corresponding to a white-noise fluctuation spectrum, is indicated 
by a dashed-dotted straight line of slope unity. The open circles and squares mark the frequencies 
1/Vg and 1/tq, characterizing the (3- and a-relaxation process, respectively. The full circles and 
full squares mark the frequencies l/t a and l/t' a , respectively. 
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FIG. 6. Evolution of the probability density P(rj,t) to find at time t the molecular axis e(t) 
with projection rj(t) onto its initial direction. The dotted lines are the initial distributions, Eq. 
(|19b| ), downscaled by a factor of 10. The oscillations around P(r],t) = are due to restricting the 
infinite sum over angular-momentum indices in Eqs. ( |l~9| ) to j < 7 (upper panel) and j < 5 (lower 
panel), respectively. 

FIG. 7. The full lines exhibit the reorientational correlators rescaled 

to c^(t) = [C^\t) — fj]/hj for two distance parameters e and the angular- momentum indices 
j = 1 — 4. The dashed lines show the /3-correlator G(t) = c a g-(t/t a ) of the hard-sphere system, 
obtained from Eqs. @, (||), and (|B"a|). 

FIG. 8. Reorientational correlators C J (t) for j = 1,2 and £ = 0.80 and £ = 0.33 for various 
distance parameters e = — 10~ x , presented as functions of log w (t / t^) . The a-relaxation-time scale 
T 3 a is defined by C^\ri) = fj/2- The horizontal lines indicate the plateaus /?. 

FIG. 9. Dipole correlator C^\t) of the dumbbell with small elongation £ = 0.33 and distance 
parameters e = — 10 _x for x = 2-5 as functions of the logarithm of the rescaled time t = t/t' a (light 
full lines). Here t' a is the second critical time scale, Eq. (|Bb|). The heavy full line is the a-relaxation 
master function C^'(t). The dotted lines show the leading-order a-scaling result plus the leading 
correction term according to Eq. ( |31b|) . The inset exhibits in a double-logarithmic plot t' a (full 
squares) and the ad-hoc scaling time r\ (open squares) for x = 1-5. The full straight line with 
slope 7 = 2.46 exhibits the power-law formula for the hard-sphere system, Eq. (|6b|). The dashed 
line interpolates the open squares for x = 1,2,3 with an effective power law exponent 7' = 1.65, 
while the dotted line is the asymptotic small-e result for t\. 
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FIG. 10. Susceptibility spectrum \" of PC as measured by incoherent neutron scattering [38] 
for q = 1.3A 1 and T = 285K (circles). The solid lines exhibit the neutron-scattering response 
of the discussed MCT model, and the dashed lines are the mere center-of-mass contributions for 
packing fractions corresponding to x = 1 and 2. The computational wave vector is q = 7.4. As 
done in Fig. |], a rescaling of the theoretical frequencies by a factor of 10 was chosen to match the 
scale of the experiment. The normalized theoretical spectra have been rescaled by a factor of 1.1. 
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